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Abstract 

The (l+l)-dimensional nonlinear boundary value problem, modeling the process of 
melting and evaporation of metals, is studied by means of the classical Lie symmetry 
| method. All possible Lie operators of the nonlinear heat equation, which allow us to 

■ reduce the problem to the boundary value problem for the system of ordinary differential 

equations, are found. The forms of heat conductivity coefficients are established when the 
given problem can be analytically solved in an explicit form. 
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1 Introduction 

It is well-known that processes of melting and evaporation of metals in the case where their 
surface is exposed to a powerful flux of energy are described by a nonlinear boundary value 
problem of the Stefan type [IrS]. In the (l+l)-dimensional case the relevant boundary value 
problem (BVP) reads as [SHI]: 



<£: Tx\ Xl(Tl) i^) =Cl(Tl)? w- < s i( f ) < ^ < f 2 w, a) 

§: |( w) f).c 2 ,r 2 ,f, ,>& W , ( 2 ) 

x = S 1 (t):X 1 (T 1 )-^ = -^H v -q(t) i (3) 
x = S 1 (t):T 1 =T v , (4) 

X: ^ = fl,(*):VT a )§ = A 1 (T 1 )^ + fF mj (5) 

x = S 2 (t) :T 1 =T 2 = T m , (6) 
x = +oo : T 2 = T , (7) 

where T v , T m , T are the known temperatures of evaporation, melting, and solid phase of 
metal, respectively; are thermal conductivities; C&, H v , H m are specific heat values per 
unit volume; q(t) is a function presenting the energy flux being absorbed by the metal; 
are the phase division boundary coordinates to be found; are the phase division boundary 
velocities; Tk(t, x) are unknown temperature fields; and index k = 1, 2 corresponds to the liquid 
and solid phases, respectively. 

In this BVP with moving boundaries, Eqs. ([[]) and 02]) describe the heat transfer process in 
liquid and solid phases, respectively, the boundary conditions fl3]) and PJ present evaporation 
dynamics on the surface Si, and the boundary conditions (J5J and (]6]) are the well-known Stefan 
conditions on the surface S 2 dividing the liquid and solid phases. Assuming that the liquid 
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phase thickness is considerably less than the solid phase thickness, one may use the Dirichlet 
condition It should be stressed that we neglect the initial distribution of the temperature 
in the solid phase and consider the process at that stage when two phases already take place. 
This means that we start to describe the process at time t = t* > when 

71 = T 2 = T s (t), 

where Ti(t) and T s (t) are non-constant functions, which are defined by the solutions of the 
problem 

The simplest realistic case of this BVP with moving boundaries occurs under the assumption 
q(t) = const when the process has a long quasistationary phase after a short transient phase 
for t G (0, t*). It means that the unknown functions Si and S 2 are linear with respect to the 
time if t > t* and S 2 — Si = const; therefore, the BVP ([I])-© can be reduced to the problem 
for ordinary differential equations by ansatz 

T 1 = T 1 (z), T 2 = T 2 (z), z = x-vt, 

where v = ^ = > is an unknown phase division boundary velocity. It turns out that the 
BVP obtained can be exactly solved in an implicit form; moreover, the solution is expressed in 
an explicit form for a wide range of functions and Ck Note that the case of constant 

values Afc and Cj-, i.e. the fact that Eqs. (JT]) and (J2J), are linear heat equations, was considered 
in the pioneering paper [8]. 

This paper is devoted to finding new reductions of the nonlinear BVP flTj)-(J7]) to the simpler 
problems and to constructing their exact solutions. The main idea is to apply the classical 
Lie symmetry method [9HTT]. In section 2, all possible Lie operators of the nonlinear heat 
equations ([I]) and ([2]), which allow us to reduce the problem to the BVP for an ordinary 
differential equation system, are found. In section 3, the forms of the coefficients arising in 
BVP (IT])— (!7|) are established when the boundary value problems obtained in section 2 can 
be analytically solved in an explicit form and the relevant exact solutions are constructed. 
Application of the exact solution obtained in the case of linear basic equations (see Eqs. ([I]) 
and (T5]) with constant coefficients) is presented to calculate the temperature fields and phase 
division boundary coordinates for the parameters, which are typical for aluminium. Section 4 
concludes the paper. 



2 Reduction of the problem to the nonlinear BVP for 
the system of ODEs 

It can be noted that BVP (JTJ) — ([T]) can be simplified if one applies the Kirchhoff substitution 

u = J Ci(Ti)dT u v = J C 2 (T 2 )dT 2 . (8) 



Substituting (jHJ) into ([I])-© and making the relevant calculations, we arrive at the equivalent 
BVP of the form 

°" ° ^di(u)^), t>0, Si(t) <x<S 2 (t), (9) 



dt dx V dx J 
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dv 
di 



dx \ dx 
Si(t) : di(u 



t>0, x> S 2 (t), 



du dS\ 



dt 



H v -q(t), 



x 
x 

' dx 1 v ' dx dt 
x = S 2 (t) : u = u 2 , v = v 2 , 

x = +00 : v — Vq, 



dx 

Si(t) : u — u\, 

ry /,\ , / \9V , . .8U dS 2 T1 

S 2 {t) : d 2 (v)— = dxiy)— + — H n 



(10) 

(11) 
(12) 

(13) 

(14) 
(15) 



where dAu) 



MTl) J _ A 2 (T2) 
Ci(Ti)' "2^; c 2 (T 2 )- 

Now one can see that BVP (M)- ffT5]) is based on the standard nonlinear heat equations 
and ffTUj) . We want to find all possible reduction of this BVP to a nonlinear BVP based on 
ODEs (not PDEs !) using the known Lie symmetry operators of the nonlinear heat equation 
(NHE). We start from theorem 1 that gives strong restrictions on the form of Lie symmetry 
operators. 



Theorem 1 A Lie symmetry operator of NHE 



du 
~di 



dx \ dx 



d{u) 7^ const 



reduces this equation together with the moving boundary conditions 

Si(t) : u — Ui, x = S 2 (t) : u 



x 



u 2 , 



(16) 



(17) 



where Sk{t) are unknown non-constant functions while u 2 > u\ are the given constants, to an 
ODE with the relevant boundary condition iff the operator in question up to local transformations 
x — )• x + xq, t — > t + to (xo € R, to ^ R) is equivalent either to 



or to 



Xi = d t + fid x , 11 G 
X 2 = 2td t + xd x . 



(19) 



Proof. The group classification of NHE f|T6|) is well-known [9]. If d(u) is an arbitrary function, 
then the maximal algebra of invariance (MAI) is generated by the basic operators (dt, d x , 2td t + 
xd x ). There are three special cases of extension of this three-dimensional algebra (see table 1). 

Table 1. Lie algebras of NHE (16). 



No. 


The form of NHE 


MAI 


1. 


u t = (e u u x ) x 




(d t , d x , 2td t + xd x , xd x + 2d u ) 


2. 


u t = (u k u x ) x ,k^ 0, • 


4 

3 


(d t , d x , 2td t + xd x , kxd x + 2ud u ) 


3. 


u t = (u~~5u x ) x 




(d t , d x , 2td t + xd x , -\xd x + 2ud u , x 2 d x - 3xud u ) 



Each NHE that admits four- or five-dimensional Lie algebra is reduced to one of those from 
table 1 by the equivalence transformations 



t = e t + t , x = e\x + x , u = e 2 u + u , 



(20) 
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where eo, to, e\, Xo, e 2 , and u are arbitrary group parameters. 

According to the Lie theory, each linear combination of the Lie symmetry operators allow us 
to reduce the relevant NHE to an ODE; however, we also need the correctly-specified reduction 
of the moving boundary condition ( TTT|) . 

Let us consider an arbitrary function d(u) and assume d(u) 7^ e u ,u k . In this case, the most 
general form of the Lie symmetry operator is 

X 3 = (\ 1 + 2\ 3 t)d t + (\ 2 + \ 3 x)d x . (21) 

Hereinafter, A with indices are arbitrary constants. It is well-known that operator (I2ip up to 
local transformations of NHE (TL61) 



x^x + x , t ->■ t + t (xo G R, t E E) , (22) 

which is a subset of (|20|) . can be reduced either to form (fl8|) (if A3 = 0) or (fl9|) (if A3 7^ 0). 
Note that the case A3 = and A1A2 = leads to the ansatz, which contradicts the moving 
boundary condition (TTTj) . 

Moreover, transformations (1221) preserve the form of condition (TTTj) (up to new notations). 
Examination of operator (lT8"j) immediately leads to the result of paper [7] because it generates 
the plane wave ansatz 

u(t, x) — U(z), z = x — fit, (23) 

which leads only to the linear form of the function S(t) in (fl7|) . 
Examination of operator (Tl"9"j) immediately leads to the ansatz 

x 

u = u{oS), uj = —=. (24) 
Vt 



Using the second formula from (|2~4"|) . one sees that the moving boundary conditions f[T7|) take 

ce th< 

s k {t) 



the form ui = : u — u k (k — 1, 2). Since these equalities must take place for arbitrary time 



t > t* , we arrive at the conditions 



u k , (25) 



Vt 

where Uk are unknown constants. 

Thus, ansatz ( 124]) reduces problem (Tl6l) and (jTTj) to the problem 

{d{u) Ul J) u + |^ = 0, (26) 

uj = u)i : u = Mi, (27) 
w = o; 2 : u = u 2 , (28) 

if the moving boundary conditions (11 7\i take the form 

x = UkVi : u = Uk- (29) 

If NHE admits four- or five- dimensional Lie algebra, then one can reduce to the form listed 
in case 1, 2 or 3 of table 1 by the substitution u = e 2 u + uo, e 2 7^ (see f l20|) ). Simultaneously, 
conditions (ri~T|) take the form 

x = Si(t) : u = e 2 U\ + uo = U\, x = S 2 (t) : u = e 2 u 2 + u = u 2 , 
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where U\ — m 2 7^ 0. Thus, we arrive at the same problem f lT6|) and (fTT]) with bars. Hereafter, 
bars are omitted so that we need to examine only three cases from table 1. 

Let us consider the first case of table 1. Here, the most general form of the Lie symmetry 
operator is 

X 4 = (Ai + 2X 3 t)d t + (A 2 + (A 3 + X 4 )x)d x + 2X 4 d u . (30) 

Clearly, we should assume A4 7^ otherwise we arrive at the previous case (see operator OTj) ). 
Depending on the values of Aj (i = 1,4), the corresponding system of characteristic equations 

dt dx du . . 

A x + 2A 3 t ~ A 2 + (A 3 + Xa)x ~ 2X~ A 

can generate only two types of ansatze 

u = F(u) + ${x), (32) 
u = F(u) + V(t), (33) 

where F(u) is a new unknown function of the known variable u(t,x) and <&(x) and \l/(£) are 
the known functions. Substituting the first of them into (|17p and dealing similar to the case of 
operator flM]), we arrive at the conclusion that x = Sk(t) — > to = u k so that (1T7I) takes the form 

u = u k : F(u) k ) + $(S k (t)) =u k . (34) 

Since these equalities must take place for arbitrary time t > t*, we arrive at the conditions 
$(x) = $(£fc(t)) = $fc = const. So ansatz (I32p can be rewritten in the form 

u = F*(u), (35) 

where F*(u) = F(oo) + § k . On the other hand, ansatz (13"5"1) can be obtained from operator (15D1 
only under condition A 4 = but we assumed A 4 7^ 0. In the quite similar way, one proves that 
application of ansatz (133]) also leads to the requirement A4 = 0. Thus, we have shown that 
there are no new reductions of problem ffT6l) and (fTT|) in case 1 of table 1. 
In case 2 of table 1, the most general form of the Lie symmetry operator is 

X 5 = (Ai + 2X 3 t)d t + (A 2 + (A 3 + kXi)x)d x + 2X A ud u . (36) 

Depending on the values of Aj (i = 1,4), the corresponding system of characteristic equations 

dt dx du , . 

Ai + 2X 3 t ~ A 2 + (A 3 + kX^)x ~ 2X 4 u ^ ' 

can generate only two types of ansatze 

u = F(u))-$(x), (38) 
u = F(u)-V(t). (39) 



Nevertheless, these ansatze differ from ansatze fl32|) and (133)) ; one can deal with them in a quite 
similar way. Finally, one arrives at the function restrictions $(x) = $0 an d ^(t) = ^0 (where 
$0 — const and \&o — const), which immediately lead to ansatz (I35p with F*(u) = F(u>)&o or 
F*(u) = F(ui)^q. Thus, there are no new reductions of problem ( )TB)) and (fT7)) in case 2 of 
table 1. 

Examination of case 3 from table 1 is rather cumbersome; however the result is still the 
same: problem (fl6)) and (jTT)) with d{u) = u~3 can be reduced to an ODE with the relevant 
boundary condition iff the Lie symmetry operator has form (1211) . 

The proof is now completed. ■ 
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Remark 1 To prove the theorem, only the structure of the corresponding Lie ansdtze was used. 
These ansdtze are listed in an explicit form in fTB{ . 



Remark 2 One easily checks that operators ( fl^j and (Q7|) reduce problem [Tty) and (77]] to the 



linear ODE with the relevant boundary condition also in the case d{u) = const. However, this 
is rather a long routine to prove that there are no new Lie symmetry operators providing the 
same reductions because the linear heat equation admits infinity- dimensional Lie algebra. 

Theorem 2 A Lie symmetry operator of NHE reduces the nonlinear BVP /fPj)-/ T73]) to a BVP 
for two ODEs with the relevant boundary conditions iff the operator in question up to local 
transformations x — > x + Xq, t — y t + to (xq G K, to £ R) equivalent either to operator [W\l or 
to ( f7P|) and the functions Sk(t), k — 1, 2, and q(t) have the correctly specified forms 

S 1 = fit + ojx, S 2 = fit + uj 2 , q(t) = q (40) 

or 

Si = uiVi, S 2 = LU 2 Vi, q(t) = (41) 

yt 

where [i and u^k = 1,2 are to-be- determined constants and go is an arbitrary positive constant. 

Proof. The proof of this theorem is based on theorem 1. One notes that the nonlinear BVP 
(JHD-fllSD contains NHE © with the boundary conditions ((ED, (HID and NHE with the 
boundary conditions (1141) . ffTB]) so that this BVP can be reduced to a BVP for ODEs only in 
the case when the given Lie symmetry operator up to local transformations x — > x + Xq, t — > 
t + t ,xo G R, to € IR is equivalent either to operator f|T8l) or to f|T9|) . 

To complete the proof, we need to check whether these operators correctly reduce the 
boundary conditions ffTTj) and f|T3|) . In paper [7J, this has been shown for operator f|T8l) and it 
was established that the phase division lines Sk{t) = fit + u^., k = 1,2, where the constants \i 
and Uk are to-be-determined. 

The application of operator ffUJl) to BVP flUl)- f|To1) leads to the ansatz 

x 

u = u(u), v = v(u), oj = —=. (42) 

yt 



To satisfy the boundary conditions (jTzJ) and f[14l) . we obtain the phase division lines of form 
(I2SD, i.e. 

Sl{t) -uj S2{t) -u (A3) 

where u^, k = 1,2 are to-be-determined constants. The direct calculations show that ansatz 
( f4"2"j) correctly reduces the boundary conditions ffTTl) and f[T51) with the restriction (|4"3|) if addi- 
tionally the energy flux is given by the function 

q(t) = ^. (44) 

Thus, substituting formulae (14"2"|) - (14"4"|) into BVP (P)- (lToT) and making the relevant simplifi- 
cations, we arrive at the BVP for two ODEs of the form 

(di(u)*O w + = °» ( 45 ) 
{di(v)v u ) u + |« w = 0, (46) 
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u = U!-. d 1 {u)u UJ = ^jH v -qo, (47) 

u = oji : u = Mi, (48) 

0J 2 

u) = u)2:d2(v)v u = d 1 (u)u u + —H m , (49) 

w = u 2 : u = m 2 , v = v 2 , (50) 

= +oo : v = Vq. (51) 



The proof is now completed. 



3 Exact solutions of the nonlinear BVP (Dp)— (H) 

It was established in the previous section that ansatz f T42l) reduces BVP ([I])-© to the BVP 
for two ODEs (|45p - (|5ip under the corresponding restrictions. Nevertheless, BVP (j45|) - (|51|) is 
much simpler than the original problem, this BVP cannot be exactly solved in the general case 
because nonlinear ODEs (|45p and (T4"6T) are integrable only in special cases. Here, we consider 
such cases in details. 

First of the all, we introduce new variables using the well-known formulae 

U V 

U= [di(u)du, V= d 2 {v)dv, (52) 



(we assume that U and V are continuous on [u*, +oo) and [i>*, +oo), respectively). 
The local substitution ( 152]) reduces BVP (I45l)- (15T|) to the form 

U uu + < ^D 1 {U)U u = 0, (53) 
VLw + ^D 2 (V)V U = 0, (54) 

u = u l :U u = y Hv ~ 9o ' ( 55 ) 

u = Wi:U=U 1 , (56) 

uj = u} 2 :Vu, = U u + Y Hm > ( 57 ) 

UJ = U j 2 :U = U 2 , V = V 2 , (58) 

uj = +00 : V = V , (59) 

where the functions U(u),V(u) and constants oji,u 2 are found and D\{U) = j^j, D 2 (V) = 
Uk v k 
U k = J di(u) du, V k = J d 2 {v) dv, k = 0,1, 2. 

Since the basic Eqs. fl53|) and ([51)1 are still nonlinear second-order ODEs, we used the 
book [13] . which is the essential extension of the classical Kamke handbook, to specify the 
integrable cases. Taking into account that ui,D\{U) and D 2 (V) should be positive (otherwise 
one obtains non- realistic equations for the given process), only three cases were separated, 
which are listed in table 2. Now one notes that the general solutions presented in table 2 can 
be applied for analytically solving BVP (I53|) -( l59|) in nine different cases. Let us consider the 
most typical of them. 
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Table 2. Solutions of ODEs of the form (53). 



No 


ODE 


General solution 


1. 
2. 
3. 


U wu + = (a = const) 
Uuu + ytjiuUu = (a = const) 
i/aw + ^e u uU u = (a = const) 


U = C 2 + d&ert (<f) 

U = d (^erf(r) + C 2 ) , w = \ fori! + de^ 2 ) 
U = f^ + C 2 ,u; = Te-i, 
C x = ln(2/ - r) - ^ - £ r 2 



Example 1: Di(U) = a 2 , D 2 (V) = b 2 (hereafter a and b are arbitrary positive constants). 



According to case 1 of table 2 the general solutions of Eqs. (1531) and (154)1 are given by the 
formulae 



(60) 
(61) 



t/^ + C^erff^ 
a \ 2 



TV 



V = C i + C 3 ^-eri 





(t) 



Hereafter, Cj(z = 1, ... ,4) are to-be-determined constants. 

Substituting solution ( 160)1 into the boundary conditions ( 156)1 and (158)) . one finds the con- 
stants C\ and C 2 : 



c/ 2 -f/i 



Co 



^ierf (gf ) - C7 2 erf {*f ) 
/ierf ) -erf (^) ' ~ z ~ erf (*f- ) -erf (^) 

Similarly, the constants C3 and C4 are found using formulae (IBT)) . ( 158)1 and (159)1 



(62) 



Vierf(^f) - 1 



C'i 



Vberf(^f) - Tj 
erf - 1 



(63) 



So, substituting formulae (162)1 and (163)) into (160)1 and (ISTl) . respectively, we find the unknown 
functions in the explicit form 



U 
V 



Ui erf (gf) - t/ 2 erf (jf) + (U 2 - erf (f) 

erf (2f) - erf (s^) 
V r erf(^f) -V2 + (V2-Vb)erf(^) 
erf(^)-l 



(64) 
(65) 



However, we also need to specify the parameters o>i and cj 2 , which allow us to find the moving 
boundaries. This can be done by substituting (IM)) and (165]) into the boundary conditions 
( 153)) and (157)) and taking into account the equality ^(erf(x)) = ^= e _:c . After the corresponding 
calculations, we arrive at the transcendent equation system 



a 



U 2 - Ut _«2A 



v^erf (*f) - erf (?f) 
b V 2 -V Q 



e ~ 



V^erf(^f) - 1 



e ^ 



H v — qo, 
U 2 -Ux 



a ^2 u 2 

neri (sf- ) -erf (^f) 6 + T m 



(66) 
(67) 



to find the parameters u\ and lj 2 - Thus, formulae ( I64p and fl65l) present the exact solution of 
BVP (J53D-059D with Dx{U) = a 2 ,D 2 (V) = b 2 . 

Here, we present the application of formulae (164]) and fl65|) for solving this BVP with the 
coefficients, which are typical for aluminium [B] . The system ( IBT)]) and ([ST]) was solved by means 
of the program MATHEMATICA 5.2: ui « 0.0127, u 2 » 0.0202. With the values wi and w 2 , 
the temperature fields for liquid and solid phases were plotted using the program MAPLE 12 
(see Fig.l). 




Figure 1: Exact solution of the problem ([I])-([7j) with parameters that are typical for aluminium: 

\, = A 2 = 240^,d = 2,7-10 6 -^,C 2 = 2,74 -10 6 -^,^ = 2, 69 ■ 10 10 H m = 

0, 17 • 10 10 ^3,T = 300 K, T m = 933 K and T v = 2793 K ; the energy flux was set q = w 
Example 2: D X {U) = D 2 (V) = 

According to case 2 of table 2, the general solutions of Eqs. (153]) and fl54]) are 

U = d erf(r) + C 2 ^j , w = ~ (2rt7 + de"^) , (68) 

F = C 3 ^erf(z/) + C 4 ^, a; = i(2z/^ + C 3 e^ 2 ). (69) 



This is important to note that the second formula in (1681) gives one-to-one correspondence 
between u; and r, because the differentiable function u(t) is strictly monotonic: 

^ = - {w + 2r^ + C 1 (-2r)e- 2 ') = * 17 > 0. (70) 
dr a \ ar J a 

Hence, the function u{r) is reversible and an inverse strictly monotonic function r = t(u;) exists 
for all u) > 0. Analogously, we prove the existence of an inverse strictly monotonic function 
v = v{uf) for the function arising in the second formula of (169]) . With the monotonic 
differentiable functions t(u>) and v(u), we transform the boundary conditions (155]) - (159]) with 
D 1 (U) = £,D 2 (V) = £to the form 

dU ( du\ l ui 

T = Ti: &\&) =T Hv - gD > (71) 

r = r l :U = U l , (72) 
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dV (du\ 1 dU (du\ 1 u 2 
r = r2) z/ = u 2 : — — = — I — I + —H m , (73) 



dv \dv ) dr \dr J 2 
r = t 2 , v = v 2 : [/ = U 2 , V = V 2 , (74) 
v = +oo : V = V Q . (75) 

Substituting solutions (jSSl) and into the boundary conditions (1721 . (J74"|) and (J75|) . one finds 
the constants Cj(i = 1, 4): 

2 £W7i y^F ^erf^-t/aerfCn) 
1 v^erf^-erfCn)' 2 2 [7 2 - U t ' 1 J 

2 y 2 -y _ y^ Voerf^)-^ 

y/n erf (z/ 2 ) - 1 2 V 2 - V 

So, substituting formulae (176|) and ( ITT)) into (|68|) and (1691) . respectively, we find the unknown 
functions in the explicit form 

U = U 1+ (Oi -ft), »=?frt/ + -L ■ e^) , (78) 

erf(r 2 ) - erf(ri) a \ v /7rerf(r 2 ) — erf(ri) / 

erffr) =JESp) (vs _ vj), w 2 („„ + e - , (79) 

erf (u 2 ) - 1 V V 71 " err (^2) - 1 / 

where r l5 r 2 , v 2 are to-be-determined parameters. To find these positive parameters, one needs 
to substitute the functions U and V into (ITT)) and (I75|) , respectively. Moreover, the condition 
(see formulae (158)) ) 

w(r 2 ) = cj(i/ 2 ) = w 2 (80) 
should be taken into account. Finally, we arrive at the transcendent equation system 

H v ) __ c ,_ \ __ c ,_ \ = U x H v ri - aq , 



V?7i 7 erf(r 2 ) - erf(n) 



F 2 erf(z/ 2 ) — 1 v^ 7 V^2 / erf(r 2 ) — erf (n 
. , 1 U 2 - Ux J\ X f 1 V 2 -V 



a \ ^/7rerf(r 2 ) — erf(n) / 6 \ y/n evt{v 2 ) — 1 

for finding the parameters T\ , r 2 and v 2 . 

Thus, BVP (J53D ([59D with D X {U) = ffr, D 2 {V) = £ has the exact solution (|7S|t ? (f7^ and 

"i = - (r^ + /^l, . e-?) , (82) 
a V v/7rerf(r 2 ) - erf(n) / 

W2 = H ( T2 c/ 2 + 4 C/ '- C/ - e-l) , (83) 
a\ y/n erf(r 2 ) - erf(n) / 

where the parameters Ti , r 2 and v 2 are found from system (|8ip . 

Example 3: Di(U) = a 2 ,D 2 (V) = b 2 e v , i.e., we consider the case when the basic equations 
( 15 3p and ( I54p are essentially different. 
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According to cases 1 and 3 of table 2, the general solutions of Eqs. (153]) and fl54j) are given by 
the formulae 



C/ = C2 + Cl ^erf(f), (84) 



V 



V = C 4 + [ u = ue-T, C 3 = \n(2g - v) - — ^— - (85) 

J g{u) 2g-u 4 

c 

where c is an arbitrary constant. 

In a quite similar way as was done in Example 2, one shows that the function is 
reversible and an inverse strictly monotone function v = v{oo) exists for all to > 0. So, using the 
monotonic differentiable function z/(w), we transform the boundary conditions fl55l) - fl59l) with 
D X (U) = a 2 , D 2 (V) = b 2 e v to the form 

lo = oji : U u = y Hv ~ 9o ' ( 86 ) 
u = u 1 :U = U u (87) 

, = „ 2: ^f^V 1 =f /„ + ^ m , (88) 



dv \dv 

u = v 2 :U = U 2 , V = V 2 , (89) 
u = + oo:V = V , (90) 

where u 2 = v 2 e~^ . 

The constants C\ , C 2 and C 4 are defined by substitution (1811) and (18"oT) into the boundary 
conditions flH7|) and f )89|) . So, we obtain after the corresponding calculations 

n a U 2 -U 1 



erf(*f ) -erf(^f)' 



= frerf(^)-£fterf(^) 
2 erf (?f) - erf (^) ' 1 J 



C 4 = ^ 2 



^2 

dv 



The constant C3 can be found using the third equation of (|85|) with z/ = z^ 2 : 

C 3 = ln(2^ 2 ) - u 2 ) - b ^u 2 2 . (92) 

Substituting (1911) into fl84l) and (1851) . we obtain 

= gt erf (flf) - erf (gj + (£/ 2 - gj erf (f) 

erf (^f ) -erf {^f) ' 1 J 

V 

/d 
w = i/ e -T. (94) 

V2 
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Finally, substituting (193 ]) . f|94|) and ( 192] into the boundary conditions (j86j) . (j88j) and fj90j) and 
making the corresponding simplification, we arrive at the transcendent equation system 

a U 2 — Ui ° 2 "i coi 

H v — go , 



0Ferf(^f ) -erf(^) 2 

oo 

dv 



V = V 2 + 

U2 



giyY 



ln(2#(z/) - v) - — — -v = \n(2g 2 - v 2 ) - -v 2 , 



2g(u) - v 4 v ^ " 2g 2 - v 2 I 

where g 2 = g(v 2 ) = f + (f F m + (f - g ) 



4 Conclusions 

In this paper, the (l+l)-dimensional nonlinear boundary value problem (1§j)-( ll5p . modeling the 
process of melting and evaporation of metals, is studied by means of the classical Lie symmetry 
method. Theorem 2 that gives all possible Lie operators, which allow us to reduce the problem 
to the BVP for the ODE system, was proved. The forms of heat conductivity coefficients 
are established when the given problem can be analytically solved in an explicit form and the 
relevant exact solutions are constructed (see the formulae in examples 1-3). 

We found that the case of a free boundary, which moves proportionally to y/i, was earlier 
established for the classical Stefan problem with one moving boundary (solidification process) 
[HJ[T5]. m the particular case, one notes that formulae ( 16%|) -(|671) with lo\ — produce the 
corresponding solution obtained in [T3] for the Stefan problem with one moving boundary. 
Similarly, formulae (I78p . (I79p and (152] are generalizations of those from [T2](see theorem 6) to 
the case of BVP with two moving boundaries. It should be noted that the authors of [T4][T5] 
did not use any Lie symmetries but an assumption that q(t) = In paper [16], a one-phase 
Stefan problem based on the linear heat equation was analytically solved using the above- 
mentioned assumption on the free boundary describing the movement of the shoreline. In [1], 
exact solutions are presented for several Stefan problems with different types of boundaries but 
with the linear basic equations. 

To the best of our knowledge, there are only a few papers devoted to constructions of exact 
solutions of nonlinear boundary value problems of the Stefan type by means of the Lie symmetry 
method. Probably, paper [17] can be quoted as the first in this direction. Nevertheless, a Stefan 
type problem seems to be a more complicated object than the standard BVP with the fixed 
boundaries; one can note that the Lie symmetry method should be more applicable just for 
solving problems with moving boundaries. In fact, the structure of such boundaries may depend 
on invariant variable (s) and this gives a possibility of reducing the given BVP to one of lower 
dimensionality. The work is in progress to construct exact solutions of a a multidimensional 
BVP using this approach. 

References 

[1] Ready J F 1971 Effects of High-power Laser Radiation (New York: Academic Press) 
[2] Rubinstein L I 1971 The Stefan problem (American Mathematical Society: Providence) 



12 



[3] Crank J 1984 Free and moving boundary problems (Clarendon Press: Oxford) 

[4] Alexiades V and Solomon A 1993 Mathematical modeling of melting and freezing processes 
(Washington: Hemisphere publishing corporation) 

[5] Cherniha R M and Odnorozhenko I G 1990 Exact solutions of a nonlinear boundary value 
problem of melting and evaporation of metals under the action of high energy flux Dopovidi 
Akad. Nauk Ukrainy (Reports of Acad. Sci. of Ukraine), ser.A 12 44-47 (in Ukrainian, 
Summary in English) 

[6] Cherniga R M and Odnorozhenko I G 1991 Studies of the processes of melting and evapo- 
ration of metals under the action of laser radiation pulses Promyshlennaya Teplotekhnika 
(Industrial Heat Technic) 13 51-59 (in Russian, Summary in English) 

[7] Cherniha R M and Cherniha N D 1993 Exact solutions of a class of nonlinear boundary 
value problems with moving boundaries J. Phys. A.: Math. Gen. 26 L935-L940 

[8] Dulnyov G N and Yaryshev N A 1967 Estimation of process of heat-mass transfer by 
interaction of the energy flux with a substance Teplofizika vysokikh temperatur (Heatphysics 
of High Temperature) 5 322-328 (in Russian) 

[9] Ovsiannikov L V 1980 The Group Analysis of Differential Equations (New York: Academic 
Press) 

[10] Olver P 1986 Applications of Lie Groups to Differential Equations (Berlin: Springer) 

[11] Bluman G W and Kumei S 1989 Symmetries and Differential Equations (Berlin: Springer) 

[12] Fushchych W I, Serov M I and Amerov T K 1993 Nonlocal ansatze and solutions of 
nonlinear system of heat equatins Ukrainian Math. J. 45 293-302 

[13] Polyanin A D and Zaitsev V F 2003 Handbook of exact solutions for ordinary differential 
equations (Boca Raton, FL: CRC Press Company) 

[14] Tarzia D A 1981 An inequality for the coefficient a of the free boundary s(t) = 2ay/i of 
the Neumann solution for the two-phase Stefan problem Quart. Appl. Math. 39 491-497 

[15] Briozzo A C and Tarzia D A 2002 An explicit solution for an instantaneous two-phase 
Stefan problem with nonlinear thermal coefficients IMA Journal of Applied Mathematics 
67 249-261 

[16] Voller V R, Swenson J B and Paola C 2004 An analytical solution for a Stefan problem 
with variable latent heat Int. J. Heat Mass Transfer 47 5387-5390 

[17] Bluman G 1974 Application of the general similarity solution of the heat equation to 
boundary- value problems Quart. Appl. Math. 31 403-415 



13 



